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A.  INTRODUCTION 


The  purpose  of  this  research  project  is  to  investigate  thermoacoustic  heat  transport 

phenomena.  The  premise  of  the  majority  of  this  research  can  be  summarized  as  follows.  While 

early  measurements  of  thermoacoustic  phenomena  qualitatively  agree  quite  well  with  the  theory 

developed  by  Wheatley1  and  Swift,  there  have  been  few  thorough,  quantitative  comparisons  of 

measurement  and  theory.  The  quantitative  studies  that  have  been  performed  previous  to  our 

3 

involvement  indicate  that  the  theory  does  not  adequately  explain  the  results. 

Thermoacoustic  phenomena  and  the  associated  devices  can  be  described  in  terms  of  two 
basic  modes  of  operation,  heat  pumps  and  prime  movers.  In  thermoacoustic  heat  pumps,  an 
acoustic  field  stimulates  the  transport  of  entropy  over  the  surface  of  an  object,  resulting  in  a  thermal 
gradient  being  established  across  that  object.  In  prime  movers,  a  thermal  gradient  is  established 
across  an  object  which  is  part  of,  or  housed  in,  an  acoustic  resonator.  At  sufficiently  high 
gradients,  sound  is  spontaneously  generated  in  the  resonator.  The  stored  thermal  energy  is 
converted  into  acoustic  energy;  the  work  output  of  the  device  is  equal  to  the  energy  in  the  acoustic 
field.  Our  study  of  thermoacoustic  processes  involves  investigating  both  of  these  fundamental 
classes  of  phenomena.  Our  previous  research  has  indicated  that  when  the  acoustic  amplitude  is 
low,  measurements  and  theory  agree  reasonably  well.  However,  when  the  amplitude  is  large, 

.there  is  considerable  disagreement.4'8  In  other  words,  in  the  linear  acoustic  regime, 

|  ^ 

|  therrhoacoustics  seems  to  be  relatively  well  understood,  whereas  in  the  nonlinear  regime,  it  is  not. 

fThe  question  now  arises  "Is  it  the  thermoacoustics  which  in  not  well  understood,  or  the  nonlinear 

acoustics?" 

In  FY  1990,  we  extended  our  previous  measurements  and  began  to  investigate  ways  to 
separate  apparent  inadequacies  of  the  theory  of  thermoacoustics  from  phenomena  caused  by 
nonlinear  acoustic  affects.  In  addition,  we  have  begun  to  consider  the  performance  of 
thermoacoustic  heat  engines  that  operate  with  an  acoustic  sound  field  that  is  not  a  perfect  standing 
wave.  Pertinent  previous  results  of  our  investigations,  along  with  a  brief  discussion  of  the 
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progress  made  during  FY  90,  follow. 


B.  SUMMARY  OF  PROGRAM 

Our  research  program  consists  of  four  projects.  The  first  project  involves  measurement  of 
thermoacoustic  phenomena  using  thermoacoustic  couples  (TACs).  The  purpose  of  the  second 
project  is  to  examine  means  of  visualizing  thermoacoustic  processes.  The  third  project  is  an 
investigation  of  heat  driven  prime  movers.  The  fourth  project  considers  the  performance  of 
thermoacoustic  heat  engines  that  operate  with  an  acoustic  sound  field  that  is  not  a  perfect  standing 
wave.  Summaries  of  progress  made  in  these  projects  during  FY  1990  follow. 

1.  TAC  Measurements 

We  have  devoted  a  large  effort  to  a  thorough  investigation  of  basic  thermoacoustic 
phenomena  using  the  simplest  class  of  thermoacoustic  engine  -  a  thermoacoustic  couple,  or  TAC. 
A  typical  TAC  is  composed  of  a  stack  of  short,  poorly  thermally  conducting  plates,  spaced  by  at 
least  a  few  thermal  penetration  depths.  The  advantage  of  using  a  TAC  is  that  its  design  allows 
assumptions  to  be  made  which  reduce  the  theory  of  thermoacoustic  heat  transport  to  its  simplest 
form.  Our  previous  research4'8  has  shown  that  the  basic  theory  agrees  relatively  well  with  the 
result  of  measurements  made  at  low  acoustic  intensities.  However,  there  is  serious  disagreement 
between  theory  and  experiment  at  higher  acoustic  intensities.  A  major  effort  prior  to  FY  90  has 
been  to  develop  a  computer  controlled  apparatus  which  would  facilitate  data  acquisition.  An 
extensive  set  of  measurements  has  been  made  with  this  apparatus,  which  allows  us  to  vary  the 
acoustic  pressure  amplitude,  the  position  of  the  TAC,  the  acoustic  frequency,  and  the  mean  gas 
pressure  and  temperature.  These  measurements  represent  one  of  the  most  complete  studies  of 
thermoacoustic  phenomena  to  date. 
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The  first  part  of  FY  90  was  devoted  to  analyzing  and  communicating  the  results  of  the  initial, 
exploratory,  phase  of  the  TAC  measurements.9  Several  conclusion  may  be  drawn  from  this  work. 
The  first  conclusion  is  that  the  thermoacoustic  effect  is  relatively,  yet  not  completely,  well 
understood  at  low  values  of  the  drive  ratio  (the  ratio  of  the  peak  acoustic  pressure  amplitude  at  a 
pressure  antinode  to  the  mean  pressure  of  the  gas)  for  the  plate  configurations  used  in  these 
measurements.  The  second  conclusion  is  that  there  are  two  distinct  regions  of  behavior  at  higher 
drive  ratios.  One  region,  which  occurs  for  drive  ratios  less  that  approximately  1.1%,  is 
characterized  by  a  linear  decrease  in  the  agreement  between  theory  and  experiment.  (See  Reference 
9  for  the  meaning  of  "linear  decrease  in  the  agreement.")  The  second  region,  which  occurs  for 
drive  ratios  greater  than  approximately  1.1%,  is  characterized  by  the  onset  of  irregularities  in  the 
temperature  difference  data  series.  Also,  the  linear  dependence  of  the  disagreement  disappears. 
Further,  because  the  agreement  is,  generally,  best  in  the  vicinity  of  the  pressure  antinode  (velocity 
node),  some  velocity  dependent  effect  may  be  the  cause  of  the  poor  agreement  at  higher  drive 
ratios,  where  the  acoustic  particle  displacement  amplitude  can  be  as  much  as  a  few  millimeters. 
There  is  also  evidence  of  an  effect  which  depends  on  the  produce  of  acoustic  pressure  and  velocity. 

In  light  of  these  results,  we  have  initiated  the  second  phase  of  the  TAC  measurements, 
intended  to  further  probe  the  discrepancies  and  behaviors  uncovered  during  the  first  phase.  The 
first  series  of  measurements  is  designed  to  investigate  whether  or  not  edge  effects  become 
important  at  the  higher  drive  ratios  and  are  the  cause  of  the  irregularities  discussed  above.  In  the 
previous  measurements,  we  measured  the  temperature  difference  developed  across  the  TAC  with  a 
series  of  thermocouples  (called  a  thermopile)  placed  along  the  front  and  back  edges  of  the  TAC,  as 
illustrated  in  Fig.  1.  Because  of  its  proximity  to  the  edge,  the  thermopile  may  be  sensitive  to 
effects  which,  though  perhaps  causing  local  deviations  in  the  temperature  profile,  do  not  affect  the 
temperature  profile  in  interior  regions  of  the  plate.  In  other  words,  although  we  observe 
considerable  deviations  between  theory  and  measurements  made  at  the  plate  edges,  the  interior  (and 
majority)  of  the  plate  may  behave  in  accordance  with  predictions. 

It  is  useful  to  examine  some  data  taken  with  a  TAC  having  an  edge-to-edge  thermopile.  An 
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Figure  1: 


Figure  2: 


The  general  design  of  a  thermoacoustic  couple  (TAC)  used  in  the  measurements.  The 
central  plate  is  instrumented  with  a  thermopile  used  to  measure  the  temperature 
difference  developed  across  the  TAC.  Four  guard  plates  surround  the  central  plate. 


0.5  1  1.5  2  2.5  3  3.5  4  4.5 

kx 

Graph  showing  the  measured  (dashed  curve)  and  the  theoretical  value  (solid  curve)  of 
the  temperature  difference  developed  across  a  TAC  as  functions  of  kx  for  helium  and  a 
drive  ratio  of  1.99%.  The  mean  gas  pressure  is  114.1  kPa.  The  acoustic  pressure 
amplitude  is  2.27  kPa  and  the  frequency  is  696  Hz.  The  theoretical  temperature 
difference  is  computed  using  Eq.  (1)  of  Ref.  9. 
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example  is  shown  in  Fig.  2,  which  illustrates  the  comparison  between  the  measured  temperature 
difference  across  a  TAC  (dashed  curve)  and  the  theoretical  value  (solid  curve).  (Refer  to  Ref.  9  for 
full  details  of  the  TAC  construction  and  experimental  parameters.)  The  TAC  is  in  helium  at  298.4 
K.  The  mean  gas  pressure  and  temperature  are  1 14.1  kPa  and  the  acoustic  pressure  amplitude  is 
2.27  kPa,  giving  a  drive  ratio  of  1.99%.  The  frequency  is  696  Hz.  The  theoretical  temperature 
difference  is  computed  using  Eq.  (1)  of  Ref.  9.  The  temperature  differences  are  plotted  as 
functions  of  kx,  where  k  is  the  wave  number  and  x  is  the  location  of  the  center  of  the  TAC.  The 
positions  of  the  particle  velocity  and  pressure  antinodes  are  indicated  on  the  graph.  Four  areas  for 
comparisons  are  obvious.  First,  the  agreement  is  quite  good  in  the  vicinity  of  the  pressure 
antinode.  The  slopes  of  the  measured  and  predicted  temperature  differences  are  approximately 
equal.  Secondly,  the  agreement  in  the  vicinity  of  the  velocity  antinode  is  poor.  The  magnitude  of 
the  slope  of  the  measured  temperature  difference  is  much  less  than  that  of  the  predicted  result. 
Third,  the  data  series  of  the  measured  temperature  difference  is  quite  irregular  and  jagged  in  the 
vicinity  of  th**  max'™""!  theoretical  temperature  difference.  Finally,  the  measured  maximum 
temperature  difference  is  less  than  and  occurs  at  a  different  position  than  the  predicted  maximum. 
In  short,  the  agreement  between  theory  and  experiment  is  poor  except  near  the  pressure  antinode. 

Figure  3  illustrates  the  theoretical  temperature  difference  as  a  function  of  kx  for  drive  ratios 
from  0.17  to  1.99%  in  helium  at  a  mean  pressure  of  1 13.8  ±  0.2  kPa,  temperature  298.4  ±0.1  K, 
and  frequency  702.8  ±  4.5  Hz.  (The  individual  values  of  the  drive  ratio,  mean  pressure, 
temperature  and  frequency  have  been  omitted,  in  order  to  reduce  clutter  on  the  graph.  The  values 
given  above  represent  the  average  and  standard  deviation  for  the  entire  data  set  The  data  presented 
in  Fig.  2  is  a  subset  of  these  data.)  This  graph  clearly  indicates  the  progression  from  a  sinusoid  to 
a  sawtooth  curve  and  the  displacement  of  the  maximum  temperature  difference  toward  the  pressure 
antinode  as  described  by  Wheatley  £t  {d.1  The  measured  temperature  differences  are  shown  in  Fig. 
4  for  the  same  experimental  conditions.  The  measurements  do  confirm  qualitatively  the  sinusoid- 
to-sawtooth  progression.  One  of  the  curves  is  dashed  to  demarcate  two  regions  of  behavior.  The 
curves  between  the  kx  axis  and  the  dashed  curve,  corresponding  to  lower  drive  ratios,  are  smooth 
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Figure  3:  Graph  showing  the  theoretical  temperature  difference  as  a  function  of  kx  for  drive  ratios 
from  0.17  to  1.99%  for  a  TAC  in  helium  at  approximately  the  same  mean  pressure  and 
frequencies  as  those  in  Fig.2.  The  data  presented  in  Fig.  2  is  a  subset  of  these  data. 


Figure  4:  Graph  showing  the  measured  temperature  difference  as  a  function  of  kx  for  drive  ratios 
from  0.17  to  1.99%  for  a  TAC  in  helium  at  approximately  the  same  mean  pressure  and 
frequencies  as  those  in  Fig.2.  The  data  presented  in  Fig.  2  is  a  subset  of  these  data. 
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and  regular.  The  curves  on  the  other  side  of  the  dashed  curve,  corresponding  to  higher  drive 
ratios,  become  jagged  and  irregular.  This  sudden  transition  occurs  in  all  of  the  data  sets  we  have 
taken  with  TACs  having  a  single  thermopile  spanning  the  TAC  edge-to-edge.  The  dashed  curve 
corresponds  to  a  drive  ratio  of  1.03%. 

In  order  to  investigate  whether  edge  effects  are  the  cause  of  any  of  these  discrepancies,  we 
have  constructed  a  TAC  with  two  thermopiles,  whose  junctions  do  not  lie  along  the  edge,  and 
repeated  some  of  the  previous  measurements.  Preliminary  results  are  shown  in  Figs.  5  and  6, 
which  show  the  temperature  difference  measured  by  the  two  thermopiles  as  a  function  of  kx.  The 
TAC  is  1.495  cm  long.  The  two  thermopiles  are  1.373  cm  and  0.682  cm  long  and  centered  on  the 
TAC.  Therefore,  the  thermocouple  junctions  lie  approximately  0.061  cm  and  0.407  cm  from  the 
edges,  respectively.  Figure  5  shows  the  data  from  the  longer  thermopile,  Fig.  6  from  the  shorter. 
The  TAC  is  in  helium  at  a  mean  gas  pressure  and  temperature  of  200  ±  1  kPa  and  295.9  ±0.1  K. 
The  frequency  range  for  the  data  is  662.8  ±  4.5  Hz.  The  drive  ratio  range  is  0.20%  -  2.70%.  The 
dashed  curve  in  both  figures  corresponds  to  a  drive  ratio  0.97%.  (Recall  that  the  dashed  curve  in 
Fig.  4  corresponds  to  a  drive  ratio  of  1.03%.)  Comparing  these  results  with  those  shown  in  Fig. 
4,  they  are  strikingly  similar.  (No  direct  comparison  is  possible  since  the  experimental  conditions 
are  different.)  Moreover,  the  transition  to  an  irregular  behavior  starts  at  a  1%  drive  ratio,  as  it  did 
with  the  edge-to-edge  design.  These  data  show  that  the  temperature  profile  of  the  interior  of  the 
TAC  behaves  the  same  as  that  measured  at  the  edge.  Therefore,  the  irregularities  previoush 
observed  are  not  an  artifact  of  placing  the  thermopile  along  the  edge.  The  irregularities  may  still  be 
the  result  of  the  edge  (e  g.,  turbulence  generated  at  the  edge),  but  they  are  not  isolated  at  the  edge. 
They  extend  at  least  half  way  to  the  center  of  the  TAC. 

We  will  continue  this  series  of  measurements  in  FY  1991  to  further  identify  the  extent  of 
these  irregularities  and  determine  their  cause.  We  also  intend  to  extend  our  measurements  to 
geometries  other  than  TACs,  in  order  to  see  to  what  extent  the  observed  behavior  is  present  in 
more  useful  thermoacoustic  engine  stacks. 
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Figure  5:  Graph  showing  the  temperature  difference  developed  across  a  TAC  as  measured  with 
the  longer  of  the  two  thermopiles.  The  TAC  is  in  helium  at  a  mean  gas  pressure  and 
temperature  of  200  ±  1  kPa  and  295.9  ±0.1  K.  The  frequency  range  for  the  data  is 
662.8  ±  4.5  Hz.  The  drive  ratio  range  is  0.20%  -  2.70%.  Tlx  dashed  curve  in  both 
figures  corresponds  to  a  drive  ratio  0.97%. 


Figure  6:  Graph  showing  the  temperature  difference  developed  across  a  TAC  as  measured  with 
the  shorter  of  the  two  thermopiles.  The  experimental  conditions  are  the  same  as  those 
in  Fig.  5.  The  dashed  curve  in  both  figures  corresponds  to  a  drive  ratio  0.97%. 
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2.  Visualization  of  Thermoacoustics 


It  is  assumed  in  the  theory  applied  to  our  TAC  measurements  that  all  parts  of  the  plate  behave 
in  the  same  manner,  that  edge  effects  and  other  complicating  factors  do  not  exist.  Therefore,  if 
such  factors  do  occur  it  is  important  to  know  about  them.  One  way  to  investigate  whether  or  not 
these  complications  play  any  significant  role  in  our  experiments  is  to  make  measurements  such  as 
those  described  above.  However,  these  measurements  with  thermopiles  give  only  the  average 
temperature  difference  across  portions  of  the  plate.  They  do  not  give  the  actual  temperature 
distribution.  Perhaps  a  better  way  is  to  find  a  method  to  "visualize"  the  temperature  distribution. 
This  notion  involves  a  completely  different  measurement  technique. 

Two  methods  of  visualizing  thermoacoustics  are  apparent.  One  technique,  proposed  by  us 
previously,  is  to  utilize  a  thermal  imaging  camera  in  order  to  photograph,  in  the  infrared,  a  plate  or 
stack  of  plates  exposed  to  an  acoustic  standing  wave.  The  result  would  be  a  real  time  image  of  the 
evolution  of  the  temperature  gradients  in  the  plate  and  could  yield  information  about  thermoacoustic 
heat  transport  which  can  not  be  gained  using  discrete  thermocouples.  A  thermal  imaging  camera  is 
available  in  the  Electro-Optics  laboratory  in  the  NPS  Physics  Department.  However,  there  are 
several  disadvantages  associated  with  its  use.  Because  the  resolution  of  these  disadvantages  could 
prove  expensive,  we  have  decided  to  initially  investigate  a  second,  though  less  sophisticated, 
techrjque  to  determine  whether  visualization  is  fruitful  at  all. 

In  this  second  technique,  strips  of  liquid  crystal,  temperature  sensitive,  sheets  replace,  for 
instance,  the  plates  of  a  TAC.  When  these  strips  are  exposed  to  a  large  amplitude  acoustic  standing 
wave,  a  visible  color  gradient  is  setup  in  the  strip.  We  have  constructed  a  transparent,  Plexiglas 
resonator  in  which  to  conduct  the  experiment  so  that  the  strips  are  readily  visible.  The  strips  have  a 
5  °C  temperature  range  over  which  the  color  changes  from  brown  to  green  as  the  temperature 
increases.  When  a  standing  wave  is  setup,  the  color  change  occurs  quickly.  Within  a  few 
seconds,  the  edge  of  the  strip  closest  the  nearest  pressure  antinode  turns  green,  while  the  other 
edge  turns  brown.  As  time  progresses,  the  green  and  brown  regions  expand  toward  the  center  of 
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the  strip,  forming  a  uniform  temperature  gradient  across  the  strip.  Eventually,  the  mean 
temperature  of  the  strip  rises  and  the  over  J1  color  pattern  shifts  towards  green.  Finally,  the  entire 
strip  indicates  that  the  temperatures  are  out  of  the  temperature  range  of  the  strip. 

Several  important  results  arise  from  these  preliminary  experiments.  The  extreme  edges 
change  color  first,  demonstrating  that  the  transported  heat  accumulates  at  the  edges  of  the  plates 
first,  as  First  pointed  out  by  Wheadey  £1 M1  A  uniform  gradient  is  established  only  after  thermal 
conduction  acts  to  destroy  this  initial  temperature  gradient.  These  observations  are  entirely 
consistent  with  the  present  understanding  of  thermoacoustics.  A  final,  yet  not  unexpected,  result  is 
that  mean  temperature  of  the  plate  increases. 

Although,  we  have  made  only  preliminary  investigations  with  the  liquid  crystal  method  of 
visualizing  theimoacoustics,  the  method  has  already  yielded  useful  qualitative  information.  Aside 
from  these  benefits,  this  technique  is  very  useful  for  demonstrating  thermoacoustics.10 

3.  Heat  Driven  Prime  Movers 


The  largest  effort  in  FY  90  has  been  devoted  to  the  investigation  of  heat  driven  prime 
movers.  There  are  two  distinct  regions  of  operation  of  a  prime  mover  -  below  and  above  the  onset 
of  self-oscillation.  We  are  investigating  each.  In  an  effort  to  establish  a  thermoacoustics  research 
pr  ogram  at  the  University  of  Mississippi,  Professor  Henry  Bass  spent  his  Sabbatical  leave  at  NPS 
to  learn  from  our  experience.  He  has  been  actively  involved  in  the  prime  mover  research.  This 
exchange  has  not  been  unidirectional.  Dr.  Bass's  expertise  in  propagation  of  sound  in  porous 
media  has  been  invaluable  to  our  efforts. 

In  order  to  investigate  the  work  output  of  a  prime  mover  below  onset,  use  is  made  of  the  fact 
that  the  conversion  of  thermal  energy  to  acoustic  energy  affects  the  net  absorption  of  acoustic 
energy  in  the  prime  mover.  As  the  temperature  difference  across  the  prime  mover  stack  (See  Fig. 
7.)  increases  from  zero  to  its  value  at  onset,  the  net  absorption  decreases  to  zero.  The  absorption 
coefficient  can  be  determined  by  measuring  the  quality  factor  Q  of  the  prime  mover,  which  in  turn. 
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Figure  7:  Diagram  of  a  thermoacoustic  prime  mover. 


Figure  8:  Graph  of  1/Q  vs  the  temperature  difference  across  the  prime  mover  stack  Filled  with 
helium  at  a  mean  pressure  of  238  kPa.  The  curve  is  the  results  of  calculations 
discussed  in  Ref.  1 1. 


can  be  related  to  the  acoustic  power  dissipated  in  the  tube.  The  dissipated  power  can  be  calculated 
using  equations  derived  by  Swift,2  which  express  the  acoustic  power  output  of  a  thermoacoustic 
engine  as  a  function  of  the  applied  temperature  difference  across  the  stack.  The  absorption 
coefficient  can  also  be  extracted  from  the  propagation  constant,  as  is  commonly  done  in  porous 
media  calculations. 1 1 

An  example  of  data  taken  below  the  onset  of  self-oscillation  is  shown  in  Fig.  8,  which  is  a 
graph  of  1/Q,  proportional  to  the  net  absorption  coefficient,  as  a  function  of  the  temperature 
difference  across  the  prime  mover  stack.12  The  data  were  obtained  by  driving  the  prime  mover  in 
the  vicinity  of  its  fundamental  resonance  and  measuring  the  frequency  response.  It  is  important  to 
point  out  that  the  acoustic  levels  involved  are  well  within  the  linear  regime.The  curve  is  based  on 
calculations  discussed  in  Ref.  11.  Though  the  overall  agreement  is  quite  good,  there  are  some 
discrepancies.  There  is  a  tendance  to  over  predict  the  absorption  at  low  temperature  differences. 
Although  not  so  obvious  from  Fig.  8,  other  data  sets  indicate  some  difficulty  in  predicting  the 
onset  temperature  The  exact  reasons  for  these  discrepancies  are  not  entirely  clear.  What  is  clear  is 
that  certain  aspects  of  the  theory  have  not  been  adequately  addressed. 

We  have  also  studied  prime  movers  above  the  onset  of  self-oscillation.13  The  prime  mover  is 
said  to  have  reached  "onset"  when  the  temperature  difference  across  the  stack  is  sufficient  for  the 
prime  mover  to  generate  and  sustain  detectable  levels  of  sound.  The  use  of  the  word  "detectable" 
is  not  meant  to  imply  that  onset  is  a  subtle  question  of  detection  thresholds.  When  onset  is 
reached,  the  observer  (and  everyone  else  in  the  room)  knows  it. 

Once  onset  of  self-oscillation  is  reached,  the  acoustic  amplitude  in  the  tube  immediately 
assumes  a  large  value,  typically  about  1%  of  the  ambient  pressure.  The  observed  waveform  is 
noticeably  non-sinusoidal.  As  more  e  ergy  is  supplied  to  the  hot  end  of  the  stack,  the  temperature 
of  that  end  increases  only  slightly  while  the  acoustic  amplitude  in  the  tube  increases  rapidly, 
reaching  as  much  as  7%  of  ambient  pressure.  Unfortunately,  as  the  acoustic  amplitude  increases, 
an  increasing  fraction  of  the  acoustic  energy  appears  as  higher  harmonics  -  harmonic  distortion 
increases.  It  is  this  distortion  that  has  occupied  our  interest.  In  order  to  understand  this 
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phenomena,  we  must  not  only  fully  understand  harmonic  generation  from  finite  amplitude  standing 
waves  but  also  any  contributing  nonlinear  generation  from  the  thermoacoustics. 

Figures  9  and  10  show  the  waveform  and  the  spectrum  of  the  sound  generated  by  the  prime 
mover  above  onset  The  gas  is  helium  at  a  mean  pressure  of  307  kPa.  The  temperature  difference 
across  the  stack  is  325  °C,  which  is  slightly  above  onset.  The  signal  exhibits  slight  distortion, 
particularly  in  the  positive  cycle.  Figure  10  shows  that  the  difference  in  spectrum  level  between  the 
first  few  modes  is  larger  than  15  dB.  Figures  1 1  and  12  show  results  for  a  temperature  difference 
of  368  °C,  the  signal  is  distorted  sharply  in  both  positive  and  negative  half-cycles.  The  difference 
between  spectrum  levels  for  the  first  few  modes  has  decreased  further  to  less  than  6  dB  and  more 
energy  has  been  spread  to  higher  modes. 

There  can  be  little  doubt  that  the  signals  at  the  higher  frequencies  are  being  generated  by  the 
fundamental.  If  we  examine  the  frequency  of  these  spectral  peaks,  they  are  integral  multiples  of 
the  fundamental  to  within  four  significant  digits  and  have  no  measurable  line  width.  Prediction  of 
harmonic  generation  (one  method  of  describing  the  non-sinusoidal  waveform)  is  not  straight 
forward  even  with  a  good  understanding  of  dissipation  mechanisms  in  the  stack.  In  a  resonant 
tube,  the  frequencies  of  the  overtones  depend  upon  dispersion  while  harmonic  generation  results  in 
exact  multiples  of  the  fundamental.  This  means  that  sound  generated  nonlinearly  will  be  generated 
at  a  frequencies  different  from  tube  resonances. 

We  have  applied  a  treatment  of  Finite  amplitude  standing  waves  by  Coppens  and  Sanders^ 
which  allows  us  to  calculate  the  the  strengths  of  the  harmonics  relative  to  that  of  the  fundamental. 
Performing  the  calculations  for  the  situation  represented  in  Figs.  9  and  10  yields  (P2/P1)  =  -20  dB 
and  (P3/P1 )  =  -32  dB.  Referring  to  Fig.  10  shows  that  the  measured  values  are  closer  to  -16  and  - 
32  dB,  respectively,  in  reasonable  agreement  with  theory.  For  the  case  shown  in  Fig.  12  the 
calculations  predict  (P2/P1)  =  -7  dB  and  (P3/P1)  =  -8  dB  while  the  measured  values  are  -7.6  and  - 
14.5  dB,  respectively.  At  even  higher  fundamental  amplitudes,  the  results  are  less  acceptable. 
Although  the  results  of  this  analysis  are  not  extremely  accurate,  they  do  serve  as  a  motivation  to 
carry  the  treatment  further.  The  reasons  for  the  discrepancies  are  not  clear. 
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Figure  11: 
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Waveform  of  the  sound  generated  by  the  prime  mover  at  a  temperature  difference  of 
368°C. 
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Figure  12:  Spectrum  of  the  waveform  shown  in  Fig.  11. 


At  this  point,  we  are  not  in  a  position  to  determine  the  relative  contributions  from  nonlinear 
propagation  and  thermoacoustic  generation  to  the  observed  waveform  distortion.  At  the  very 
minimum,  such  a  determination  will  require  an  accurate  model  for  dispersion  in  the  prime  mover. 
This  area  is  receiving  attention  now.  Beyond  that,  the  equations  which  give  rise  to  acoustic  gain 
have  been  linearized.  A  numerical  solution  of  the  coupled  momentum  and  energy  equation  might 
identify  significant  nonlinear  terms  in  the  coefficient  of  gain. 

One  aspect  of  prime  movers  that  has  been  ignored  until  now  is  their  ability  to  go  into  self 
oscillation  for  more  than  one  mode  simultaneously.  In  our  prime  mover,  the  temperature 
difference  required  for  onset  for  the  fundamental  mode  is  approximately  320  °C,  while  that  for  the 
first  overtone  (the  second  mode)  is  465  °C.  By  introducing  a  disk  with  a  hole  of  the  right 
dimensions  in  it  into  the  resonator  at  the  right  place,  we  are  able  increase  the  onset  temperature 
d  fference  for  the  fundamental  without  affecting  that  for  the  second  mode.  This  allows  us  to  excite 
both  modes  simultaneously.  If  the  attenuator  (disk)  is  not  used  for  the  fundamental,  we  are  unable 
to  supply  enough  power  to  the  prime  mover  to  reach  onset  for  the  second  mode.  The  fundamental 
transports  too  much  heat!  As  one  might  imagine,  the  spectrum  of  the  waveforms  generated  is  quite 
rich,  containing  harmonics  of  both  modes  as  well  as  sum  and  difference  frequencies.  One  reason 
for  exciting  both  modes  is  that  it  gives  a  measure  of  the  dispersion  in  the  prime  mover,  a  quantity 
needed  to  understand  the  harmonic  generation  discussed  above.  It  also  provides  another  test  bed 
for  theory.  We  intend  to  pursue  these  measurements  further. 

4.  Jjon-oerfect  StandinE  Wave  Heat  Engines 

It  is  of  interest  to  consider  the  performance  of  thermoacoustic  heat  engines  that  operate  with 
an  acoustic  sound  field  that  is  not  a  perfect  standing  wave.  In  such  a  system  the  phase  between 
acoustic  pressure  and  acoustic  velocity  is  somewhere  between  zero  and  90  degrees.  This  is  a 
difficult  theoretical  problem  to  consider  in  any  generality,  given  the  complexity  of  the  problem. 
The  complexity  can  be  seen  in  the  large  parameter  space  required  to  describe  any  specific  engine. 
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There  are  gas  parameters:  sound  speed,  Prandtl  number,  Cp/cv,  and  mean  pressure.  There  are 
stack  geometry  parameters:  plate  separation,  plate  thickness,  and  stack  length.  There  are  operating 
conditions:  acoustic  pressure  amplitude,  velocity  amplitude,  their  relative  phase,  frequency,  and 
the  heat  load  on  a  heat  pumping  stack  or  the  work  load  on  a  prime  mover  stack.  These  parameters 
could  be  expressed  differently  and  equivalently. 

The  differential  equations  of  N.  Rott^.lfi  are  capable  of  describing  the  system  with  all  the 
detail  that  we  need.  The  only  significant  limitation  of  these  equations  is  that  they  are  linear  and  are 
expected  to  compare  well  with  experimental  results  only  at  low  acoustic  amplitudes.  In  some 
ways,  these  equations  require  that  too  many  parameters  and  initial  conditions  be  specified.  Some 
of  die  parameters  can  be  eliminated  by  realizing  that  not  all  of  them  are  independent,  and  that  the 
parameter  space  can  be  reduced  to  a  minimal  independent  set  by  using  fairly  obvious 
normalizations.  Two  examples  are:  expressing  the  distances  in  the  wave  propagation  direction  in 
terms  of  wavelength  units,  and  expressing  distances  in  the  orthogonal  diffusion  direction  in  terms 
of  penetration  depth  units.  However,  the  reduced  set  of  parameters  is  still  complex,  and  this 
interferes  with  the  study  of  generalized  questions  such  as  the  significance  of  acoustic  phasing. 

In  an  effort  to  simplify  the  equations  of  Rott,  Swift2  makes  additional  approximations.  The 
first  approximation  is  that  the  length  of  the  stack  is  infinitesimally  short,  and  the  second  is  that  the 
spacing  between  plates  in  the  stack  is  at  least  several  boundary  layers  in  size,  such  that  the 
diffusive  part  of  the  thermal  wave  is  nearly  attenuated  to  zero  at  the  midpoint  between  two  plates. 

The  approximation  of  short  stack  length  is  a  reasonable  way  of  simplifying  the  system,  at 
least  for  our  present  purposes.  For  the  purpose  of  analyzing  many  real  or  potential  experimental 
systems,  the  approximation  is  not  very  accurate.  However,  for  this  study,  there  is  not  any 
substantial  loss  of  information.  For  example,  a  finite  length  stack  could  be  analyzed  accurately  by 
modeling  it  as  many  infinitesimal  stacks  placed  in  series.  In  fact,  this  is  roughly  equivalent  to  the 
method  used  by  numerical  algorithms  for  solving  ordinary  differential  equations. 

The  boundary  layer  approximation  is  a  severe  limitation  for  this  study.  The  parameter  space 
of  interest  lies  somewhere  between  two  extremes.  One  extreme  is  the  "classic"  thermoacoustic  heat 
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engine  as  described  by  Wheatley1  and  Swift2.  This  engine  operates  in  a  sound  field  that  has 
acoustic  pressure  and  acousdc  velocity  in  quadrature,  and  the  spacing  between  plates  of  the  stack  is 
at  least  several  penetration  depths.  The  other  extreme  is  the  Stirling  cycle  engine,  which  has  the 
dynamic  pressure  and  velocity  in  phase,  and  the  spacing  between  plates  of  the  "regenerator"  much 
less  than  a  penetration  depth.  Dearly,  the  boundary  layer  approximation  is  not  suitable  for  an 
analysis  of  engines  that  approach  the  Stirling  cycle  extreme. 

Our  approach  will  be  to  use  Swift's17  equations  initially,  and  then  use  the  Rott15*16  equations 
with  the  short  stack  approximation  for  the  actual  study  of  acoustic  phasing.  Swift's  equations  for 
acoustically  stimulated  enthalpy  flow  H2  and  acoustic  work  flow  W2  are. 
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The  initial  use  of  Swift's  equations  will  assist  in  developing  the  analysis  structure  with  a  simplified 
set  of  equations  before  moving  to  the  more  complicated  Rott  equations.  Also  the  initial  analysis 
can  be  compared  to  some  of  Swift's  own  results  for  purposes  of  error  checking,  and  the  final 
analysis  can  again  be  compared  to  the  initial  analysis  in  limiting  cases. 

Another  technique  that  can  be  applied  to  our  problem  is  the  graphical  display  of  engine  data. 
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using  more  than  one  independent  parameter,  such  as  a  contour  plot  of  two  independent  parameters. 
Figures  13  and  14  are  two  such  contour  plots  based  on  the  above  equations  and  a  purely  reactive 
standing  wave  sound  field.  Figure  1 3  shows  contours  for  a  prime  mover  engine  where  the  solid 
contours  represent  efficiency  relative  to  Carnot  (perfect)  performance  and  the  dashed  contours 
represent  dimensionless  power  density  in  the  form  of  acoustic  work,  for  a  stack  having  a  fixed 
temperature  difference.  Figure  14  shows  contours  for  a  heat  pump  where  the  solid  contours 
represent  COP  relative  to  Carnot  performance  and  the  dashed  contours  represent  dimensionless 
power  density  in  the  form  of  heat  flow.  The  parameters  chosen  for  these  plots  are  typical  for 
engines  utilizing  a  10  bar  helium  working  fluid  with  modest  acoustic  amplitudes.  The  two 
independent  parameters  are  kx,  which  is  the  normalized  position  of  the  stack  in  the  standing  wave 
and  T,  which  is  the  normalized  temperature  gradient  in  the  stack. 

A  number  of  interesting  observations  can  be  made  from  the  contour  plots.  Most  engine 
designs  would  not  choose  an  operating  point  that  coincides  with  the  region  of  highest  power 
density  (since  the  efficiency  is  zero  there),  nor  would  one  choose  an  operating  point  that  coincides 
with  the  region  of  highest  efficiency  (since  the  power  density  is  quite  low  there).  Instead,  a  more 
likely  choice  of  operating  point  is  one  that  lies  somewhere  on  the  curve  of  best  compromise 
between  power  density  and  efficiency.  This  curve  (not  shown  in  the  figures)  is  the  one  that 
indicates  the  path  in  the  T  -  kx  plane  such  that  a  given  sacrifice  in  efficiency,  gains  the  most  in 
power  density  (or  vice-versa). 

The  work  that  remains,  is  to  produce  similar  contours  for  engine  performance  that  are  based 
on  the  Rott  equations  directly,  and  then  vary  the  acoustic  phasing.  Also,  numerically  producing 
the  curve  of  best  compromise  on  the  contour  plots  would  be  useful. 
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Prime  Mover  Power  ic  Efficiency  Contours  (3) 


Figure  13:  Contour  plot  for  a  prime  mover.  The  solid  contours  represent  efficiency  relative  to 
Carnot  (perfect)  performance  and  the  dashed  contours  represent  dimensionless  power 
density  in  the  form  of  acoustic  work,  for  a  stack  having  a  fixed  temperature  difference. 


Refrig  COPR  &  Cooling  Power  Contours  (4) 


Figure  14:  Contour  plot  for  a  heat  pump.  The  solid  contours  represent  COP  relative  to  Carnot 
performance  and  the  dashed  contours  represent  dimensionless  power  density  in  the 
form  of  heat  flow. 
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